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Summary 

A  new  low-speed  aerodynamic  database  for 
the  F-16/VISTA  was  developed  from  low-speed  data 
sets  obtained  from  NASA  Langley's  30x60  facility  and 
20-ft  vertical  wind  tunnel  and  NASA  Ames'  12-ft 
pressure  tunnel.  This  new  model  was  designed  to  be 
continuous  from  -80°  to  +90°  angle  of  attack,  -30°  to 
+  30°  of  sideslip,  and  to  incorporate  all  sideslip  and 
control  effectiveness  interactions  in  a  fully  non-linear 
structure.  The  aerodynamic  data  set  also  includes 
wind-axis  and  body-axis  damping  terms  and  is 
configured  to  utilize  these  data  separately  for  high- 
fidelity  modeling  of  dynamic  conditions.  The  purpose 
of  this  effort  is  to  demonstrate  the  benefits  of  a 
completely  nonlinear  aerodynamic  data  base  to 
capture  the  dependencies  that  occur  while 
maneuvering  at  high  angle  of  attack.  This  model  was 
validated  against  a  set  of  high  angle  of  attack  flight 
data  from  the  F- 16  VISTA/MATV  flight  test  program. 
Minimal  refinements  to  the  data  set  were  required  to 
produce  an  excellent  correlation  between  the 
simulation  model  and  the  flight  test  results  for 
maneuvers  in  and  above  the  stall  regime. 

Introduction 

With  extensive  backgrounds  in  modeling 
nonlinear  aerodynamic  characteristics,  Bihrle  Applied 
Research  was  tasked  by  the  Wright  Lab  Flight  Control 
Division  to  develop  a  new  F-16  aerodynamic  model 
with  the  philosophy  of  incorporating  all  known 
nonlinear  static  and  dynamic  characteristics  from 
available  wind  tunnel  data.  This  new  F-16  (and 
VISTA)  low-speed  simulation  aerodynamic  model 
was  developed  to  assess  the  impact  of  incorporating 
these  non-linear  effects  on  high-angle  of  attack 
simulation  to  analyze  the  VISTA/MATV  flight  test 
data  and  was  developed  to  be  amenable  to  database 
updates.  This  was  accomplished  by  applying  the 
philosophy  of  using  all  available  wind  tunnel  data  and 
letting  the  test  data  functionalities  and  all  interactions 
be  modeled  as  continuous  functions  in  the  simulation 
data  set.  This  low-speed  model  is  entirely  non-linear 


and  continuous  from  -80°  to  +90°  angle  of  attack  with 
non-symmetric  mechanization  of  sideslip  effects  for 
the  basic  aircraft  characteristics.  All  sideslip  and 
control  surface  effectiveness  interactions  were 
modeled  as  dictated  by  the  wind  tunnel  data.  While 
benefiting  from  hindsight  that  earlier  models  lacked,  it 
was  hoped  that  by  starting  with  a  carefully  structured 
non-linear  model  merging  the  best  available  data 
would  result  in  a  simulation  that  would  more  closely 
model  the  VISTA’s  post-stall  characteristics. 

This  paper  first  reviews  the  model  structure 
and  documents  the  aerodynamic  data  source  used  to 
assemble  the  data  set.  The  subsequent  discussion 
examines  the  methodology  used  to  correlate  and 
validate  the  model  against  flight  test  data  as  well  as 
the  rationales  for  the  model  changes.  Finally, 
simulation  runs  using  flight-recorded  stick  inputs  will 
be  used  to  demonstrate  the  low-speed  model's  fidelity 
in  the  stall  and  post  stall  region. 

Low-Speed  Aerodynamic  Model 

The  basic  static  F-16’s  low-speed 
aerodynamic  model  was  assemble  primarily  from  the 
NASA  Langley  30x60  low-speed  tests  for  the  upright 
angle  of  attack  range  and  the  NASA  Ames  12-ft 
pressure  tunnel  data  for  the  inverted  characteristics  as 
illustrated  in  Figure  1.  The  30x60  wind  tunnel  data 
range  extended  from  -5°  to  +60°  angle  of  attack  and 
out  to  ±20°  of  sideslip.  Because  of  the  known  non¬ 
linear  variations  in  the  aerodynamic  characteristics 
beyond  20°  of  sideslip,  the  low-speed  wind  tunnel  data 
from  the  Langley’s  20-ft  vertical  wind  tunnel  were 
used  to  extend  the  sideslip  functionality  to  ±30°  for  all 
six  body-axis  force  and  moment  coefficients.  In 
addition,  the  same  data  set  was  used  to  extend  the 
angle  of  attack  from  60°  to  90°. 

For  the  30x60  upright  wind  tunnel  data  set, 
control  effectiveness  runs  were  conducted  on  the  basic 
F-16  configuration  with  primarily  full  trailing  edge 
deflections.  Additionally,  since  a  comprehensive  set  of 
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control  surface  deflections  had  been  tested  at  the 
NASA  Langley’s  20-ft  wind  tunnel,  these  data  (from 
Reference  1  to  5)  were  used  to  define  the  nonlinear 
variation  in  control  authority  as  a  function  of  surface 
deflections  and  combined  with  control  effectiveness 
data  from  the  30x60  data  to  result  in  a  very  high- 
fidelity  control  effectiveness  model.  NASA  Ames’ 
inverted  wind  tunnel  data  were  used  to  define  the 
control  effectiveness  model  at  inverted  angles  of 
attack.  Both  the  Ames  inverted  and  the  Langley  data 
sets  exhibited  sideslip  dependencies,  and  the 
characteristics  were  incorporated  as  well. 

While  the  previous  aerodynamic  data  sets 
incorporate  only  linearized  body-axis  damping  * 
derivatives  to  model  the  F-16’s  dynamic 
characteristics,  the  new  low-speed  model  use  both 
body-axis  damping  derivatives  as  well  as  wind  axis 
rotary  data.  These  data  were  mechanized  separately 
and  combined  in  the  simulation  using  the  techniques 
proposed  by  Kalviste  (Reference  6).  For  the  wind- 
axis  damping  terms,  traditionally  referred  to  as  the 
rotary  balance  data,  the  extensive  test  results  from 
NASA  Langley  20-ft  Vertical  Tunnel  were  used. 
Nonlinear  variations  in  the  forces  and  moments  due  to 
rotation  rates,  sideslip  angles  and  control  surface 
deflections  were  also  mechanized  in  the  model. 

NASA  Langley’s  30x60  forced  oscillation  rig 
provided  the  data  for  the  body-axis  damping 
derivatives. 

To  account  for  the  differences  in  the  external 
features  of  the  VISTA  configuration,  incremental 
aerodynamic  effects  were  applied  to  the  baseline  F-16 
data  base  to  model  the  VISTA  aerodynamic 
characteristics.  Rather  than  incorporating  these  effects 
as  increments  to  the  baseline  F-16  data  set,  BAR 
elected  to  add  increments  to  the  appropriate 
components  in  the  model  that  resulted  in  replacement 
data  tables.  This  approach  maintained  the  same  data 
structure  in  terms  of  the  independent  variables  and 
their  break  points  between  the  baseline  F-16  and  the 
VISTA  configuration  while  significantly  speeds  up  the 
process  of  sim-to-flight  correlation.  These  increments 
are  derived  primarily  from  Wright  Lab’s  linearized 
data  which  were  obtained  from  Lockheed/Fort  Worth 
as  part  of  the  VISTA  simulation  package. 

Additionally,  large  sideslip  wind  tunnel  data  acquired 
for  the  Peace  Marble  II  configuration  at  NASA 
Langley's  20  ft  wind  tunnel  were  merged  and  utilized 
at  upright  angles  of  attack  as  illustrated  in  Figure  2. 


The  typical  for  all  axes  for  the  final  low- 
speed  F-16  model  can  be  illustrated  in  the  build-up  of 
aerodynamic  rolling  moment  shown  below: 

cijoui  -CUW) 

+  DCl*  (a,p*SGN(8A),|5A|)  x  SGN(SA) 

+  DCldr  (a, p *  SGN (SR),|8R|)  x  SGN(SR) 

+  DClddh(a,P*SGN(8dh),8h,|8dh|)  x  SGN(8dh) 

+  DCldhR(a,(nb)/(2V),8Hwindw„d) 

+  ClROT(a,(nb)/(2V)*SGN(P),  |p[)  x  SGN(P) 

+  DClp  (a)  x  (P  b)/(2V) 

mod 

+  DC1r  (a)  x  (R  b)/(2V) 

mod 

The  total  lateral  characteristics  are  composed 
of  the  basic  lateral  stability,  contributions  from  control 
surfaces,  rotary  (wind  axis  rotation),  and  forced 
oscillation  (body-axis  damping)  terms.  Known 
nonlinear  aerodynamic  characteristics  and  coupling 
effects  between  sideslip  angle  and  control  surface 
deflections  from  wind  tunnel  data  were  incorporated 
in  the  data  tables. 

Simulation  /  Flight  Test  Correlation  -  "Overdrive" 

For  the  validation  of  the  new  low-speed 
aerodynamic  model,  time  histories  of  the 
VISTA/MATV  flight  test  results  were  obtained  from 
Wright  Lab.  These  flight  data  included  test  points 
taken  from  MATV  flights  no.  89,  94,  and  100  as  well 
as  some  functional  check  flight  departures  of  the 
VISTA  configuration  with  thrust  vectoring  disabled. 

One  of  the  flight  test  validation  tools  used  in 
the  simulation  environment  used  in  this  analysis, 
"Overdrive",  allows  the  validation  of  simulation 
aerodynamic  database  against  flight-extracted  data 
using  the  process  illustrated  in  Figure  3. 

At  each  time  slice,  extraction  of  aerodynamic 
moment  coefficient  from  flight-recorded  time  history 
occurs  as  shown  on  the  right  side  of  Figure  3. 

Angular  rates  are  numerically  differentiated  to  obtain 
the  angular  acceleration  of  the  vehicle.  After  the 
removal  of  the  inertial  effects  and  other  external  forces 
and  moments  such  as  thrust  vectoring,  the  remainder 
is  non-dimensionalized  to  calculate  the  aerodynamic 
moment  coefficients  during  flight.  The  translational 
acceleration  signals  from  the  flight  time  history  were 
also  processed  at  each  time  step  to  extract  the 
aerodynamic  forces  acting  on  die  aircraft 

Also  at  each  time  step,  flight-recorded  states, 
such  as  angle  of  attack,  angle  of  sideslip,  control 
surface  positions,  etc.,  are  used  to  exercise  the 
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aerodynamic  model  in  accordance  with  the 
aerodynamic  model  specification  discussed 
previously.  Each  aerodynamic  model  element  (Le., 
pitching  moment  due  elevator,  pitching  moment  due 
to  flap,  etc.)  are  stored  and  summed  as  prescribed  in 
the  aerodynamic  model. 

By  overplotting  the  model  predicted 
coefficients  with  the  flight-extracted  total  coefficients, 
differences  can  be  easily  identified.  Correlating  the 
discrepancies  with  the  excitation  of  individual 
elements  and  parameters  from  the  flight  time  history 
helps  to  isolate  the  weaknesses  in  the  aerodynamic 
model. 

It  should  be  emphasized  that  there  is  no 
integration  during  an  overdrive  run.  The  states  are 
completely  restricted  to  the  values  from  the  time 
history.  The  advantage  of  this  approach  is  that  there 
no  propagation  of  error  over  time.  Any  differences 
between  the  model-predicted  and  the  flight-extracted 
values  are  strictly  the  result  of  local  error.  This 
"analog-matching"  approach  alleviates  the  problem 
associated  with  the  propagation  of  error  over  time;  any 
noise  or  drop-out  of  signals  would  not  affect  the 
analysis  of  the  subsequent  events  in  the  time  history. 
Additionally,  flight  control  system  is  by-passed  in  this 
methodology  thus  avoiding  any  confusion  between  the 
error  caused  by  aerodynamic  modeling  and  the  error 
from  the  flight  control  system  model.  Once  the  user  is 
satisfied  with  aerodynamic  model,  the  entire 
simulation  (including  aerodynamic  model  and  flight 
control  system  modeling)  can  be  validated  by  running 
the  simulation  in  open-loop  controlled  by  the  flight- 
recorded  control  surface  or  pilot  input  time  history. 

An  example  of  this  methodology  is  the 
overplot  of  the  sim-predicted  and  flight-extracted 
yawing  moment  coefficient  during  flight  89-02068  as 
shown  in  Figure  4.  The  test  vehicle  was  stabilized  at 
35°a  while  performing  a  steady  sideslip  to  -20  °p  as 
shown  in  the  lower-left  plot. 

The  first  (upper-left)  plot  compares  the 
simulation-predicted  total  yawing  moment  against  the 
flight-extracted  total  yawing  moment.  The  close 
correlation  between  the  two  curves  validates  the 
aerodynamic  yawing  moment  model  for  the  baseline 
VISTA  configuration  at  these  flight  conditions. 

Individual  elements  of  the  simulation’s 
yawing  moment  can  also  be  plotted  (upper-right)  to 
help  identify  the  source  of  any  mismatch  between  the 
sim-predicted  and  the  flight-extracted  aerodynamic 


yawing  moment  coefficient.  Here,  the  static 
directional  stability  (Cn_basic)  and  the  incremental 
yawing  moment  due  to  rudder  deflection  (DCn_dr)  are 
shown. 

Finally,  the  time  history  of  the  flight  test 
control  surfaces  are  shown  in  the  4th  (lower-right) 
plot. 

VISTA’s  Post-Stall  Lateral  Characteristics 

A  positive  static  offset  was  observed  in  the 
wind  tunnel  data  and  was  incorporated  in  the  initial 
aerodynamic  model.  However,  the  flight  test  data 
revealed  a  negative  static  rolling  moment  bias,  thus 
prompting  the  need  to  mirror  the  nonlinear  static 
lateral  stability  (CI_basic)  along  with  the  zero  sideslip 
offset  to  reflect  the  flight  test  vehicle’s  lateral 
characteristics  in  the  35°  to  45°a  range. 

The  flight  test  points  in  the  35-37.5°  region 
indicated  the  VISTA  configuration’s  ability  to 
stabilize  at  large  negative  sideslip  angles  (in  excess  of 
-15°p)  with  little  roll  control  input  while  consistently 
demonstrated  the  vehicle’s  inability  to  stabilize  at 
positive  sideslip  angles  greater  than  +5°p.  This 
observation  lead  to  the  reduction  of  VISTA’s  lateral 
stability  in  the  35-37.5°  angle  of  attack  range  to  nearly 
neutral  for  negative  sideslip  angles  as  shown  in  Figure 
5.  While  the  original  wind  tunnel  data  exhibited  some 
asymmetry  due  the  sign  of  sideslip  induced  by  the 
offset,  the  flight  extracted  asymmetry  was  more 
pronounced,  albeit  not  unreasonably  so  based  on  other 
wind  tunnel  results  of  the  F-16  configuration. 

Because  of  the  nearly  neutral  lateral  stability  at 
negative  sideslip,  no  rolling  moment  from  the  roll 
control  surfaces  (i.e.,  flaperons,  and  differential  tails) 
were  required  to  maintain  a  wings-level  attitude  when 
attempting  to  achieve  large  negative  sideslip 
excursions  (Figure  6).  However,  the  stable  lateral 
characteristics  at  positive  sideslip  range  coupled  with 
the  reduced  roll  authority  available  to  the  flight 
control  system  in  the  stall  region  permitted  less  than 
5°  of  stabilized  sideslip  excursion  (Figure  7). 

Post-Stall  Normal  Force  Characteristics 
The  correlation  of  normal  force  coefficient 
between  the  flight  test  data  and  the  model-predicted 
values  revealed  significant  differences  in  the  region 
between  22.5  to  50°ot.  This  difference  is  clearly  seen 
in  Figure  8  where  the  comparison  of  sim-predicted 
and  flight-extracted  normal  force  coefficient  are 
shown  in  the  upper  left  plot.  The  individual 
components  of  the  simulation  normal  force  coefficient 
are  plotted  on  the  upper-right  hand  comer  while  the 
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horizontal  stabilators  and  the  pitch  nozzle’s  deflection 
time  history  are  shown  in  the  lower  right  hand  comer. 

Besides  the  basic  normal  force,  no  other 
normal  force  components  were  significantly  excited 
and  virtually  neutral  stabilator  and  pitch  nozzle 
deflections  were  seen  during  the  flight.  The 
differences  can  only  be  attributed  to  the  basic  normal 
force  (CNf_basic)  of  the  VISTA  configuration.  Based 
on  thorough  evaluation  with  many  sets  of  time 
histories  in  flights  89,  94  and  100,  the  basic  normal 
force  coefficient  is  reduced  in  the  stall/post-stall 
region  to  model  the  values  seen  in  the  flight  test  data. 
A  comparison  of  the  VISTA’s  baseline  normal  force 
coefficient  between  the  original  and  the  revised  data 
set  is  shown  in  figure  9  for  0°p.  The  difference  of  0.2 
to  0.3  is  seen  throughout  the  stall  region  across  a  wide 
range  of  sideslip  angles.  It  should  be  noted  that  this 
discrepancy  between  flight  test  and  wind  tunnel 
normal  force  has  been  observed  and  accounted  for  in 
other  flight  test  evaluations  (e.g.  Reference  7). 

Validation  and  Verification  of  Low-speed  Model 

Once  the  aerodynamic  model  was 
satisfactorily  validated  against  the  available  flight  test 
results  using  the  overdrive  methodology,  simulations 
of  the  maneuvers  performed  during  flight  were 
conducted.  Pilot  input  time  history  recorded  during 
the  flight  test  were  used  as  inputs  to  the  flight  control 
system  model  (Reference  8)  while  exercising  the  new 
aerodynamic  model.  In  Figure  9,  the  new  model 
responses  are  plotted  with  the  results  of  the  original 
Wright  Lab  aerodynamic  model  simulation  and 
compared  against  the  recorded  behavior  of  the  flight 
test  vehicle.  Since  only  the  aerodynamic  model 
differed  between  the  two  sets  of  simulation,  the  close 
match  clearly  demonstrates  the  improved  fidelity  of 
the  new  aerodynamic  model. 

Conclusions 

This  paper  has  documented  the  method  and 
data  used  in  the  development  of  a  new  low-speed 
aerodynamic  model  for  the  F- 16/VISTA 
configuration.  The  approach  taken  was  to  use  the 
underlying  data  dependencies  as  dictated  by  the  wind 
tunnel  data  to  build  a  continuous,  non-linear  model. 
The  use  of  high  angle  of  attack  flight  test  data  as  a 
simulation  validation  and  database  update  tool  was 
demonstrated  using  an  interactive  method  to  quickly 
assess  many  flight  test  points.  The  satisfactory 
correlation  between  the  flight-extracted  aerodynamic 
coefficients  and  the  model-predicted  values  clearly 
demonstrates  that  aerodynamic  characteristics  in  the 
stall  and  post-stall  region  can  be  accurately  modeled. 


This  was  accomplished  by  the  preservation  of  non¬ 
linear  characteristics  from  the  wind  tunnel  results  in 
the  simulation  model  structure.  Further,  the  clearly 
defined  and  continuous  aerodynamic  model  structure 
permitted  expeditious  identification  of  the  components 
that  caused  any  mismatch  between  the  model  and  the 
flight  results  and  allowed  easy  model  update  as 
dictated  by  flight  test  extraction. 
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Figure  3.  "Overdrive"  block  diagram. 
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Overdrive"  output  comparing  flight-extracted  yawing  coefficient  against  sim-predicted  value  for  FLT  89-02068. 


Figure  6.  Comparison  of  sim-predicted  Cl  against  flight-extracted  values  for  FLT  89-02068 


Figure  8.  Comparison  of  sim-predicted  normal  force  coefficient  (CNf)  against  flight-extracted  values  for  FLT  89-84810 
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Figure  9.  Comparison  of  the  original  and  the  modified  normal  force  coefficient  for  the  F- 16  data 
base. 


